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Abstract 

The amplitudes of the — > 3tt and if — > 2tt decays are ex- 
pressed in terms of different combinations of one and the same set of 
CP-conserving and CP-odd parameters. Extracting the magnitudes 
of these parameters from the data on K — > 2ir decays, we estimate 
an expected CP-odd difference between the values of the slope pa- 
rameters g + and g~ of the energy distributions of "odd" pions in 
K + — > 7T + 7r + 7r~ and K~ — ► 7r~7r~7r + decays. 
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1. Introduction 



The observation of CP effects in K ± — > 3n decays would allow to understand 
better how the mechanisms of CP violation work. 

Now the Collaboration NA48/2 began a search for such effect with accu- 
racy 6((=£) < 2 x fCT 4 . 

Contrary to the case of Kl — > 2tt decay, where CP violates both in AS = 
2 and AS = 1 transitions, in the K ± — > 3ir decays, only the last (the so-called 
"direct" ) CP violation takes place. Experimentally, an existence of the direct 
CP violation in K L — > 2ir decays, predicted by Standard Model (SM) and 
characterized by the parameter e', is establised: e'/e = (1.66 ± 0.16) x 1CT 3 . 

What is expected for CP effects in K + — > 7r ± 7r ± 7r =F decay? To give an an- 
swer, it is necessary to understand the role of the electroweak penguin (EWP) 
operators in both decays and get rid of the large uncertainties usual for the 
theoretical calculations. The real scale of these uncertainties is characterized 
by the following predictions obtained before the above experimental result: 

- = (17^) x icr 4 [i], - = (1.5 -T- 31.6) x irr 4 [2] 

£ £ 

. To avoid the uncertainties arising in the theoretical calculation of the 
ingredients of the theory, we use the following procedure. We express the 
amplitudes of K — > 2n and K ± — >• 37r decays in terms of one and the same 
set of parameters, and calculating g + — g~ , we use the magnitudes of these 
parameters extracted from data on K — > 2n decays. 

2. The scheme of calculation 

AS = 1 non-leptonic decays is based on the effective Lagrangian 
L(AS = 1) = V2G F sin 6c cos 6c ^ CjOj , (1) 



Oi = s L -f^d L ■ ul^^ul - sl^^ul ■ u L ^^d L ({8/}, AI =1/2); (2) 
+ 2s Llll d L ■ s Llll s L ({84, AI = 1/2); (3) 



A theory of 
[3] 



where 
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- 3s l7/ A • s Ll ,s L ({27}, A/ = 1/2); (4) 



- s Ll ^d L ■ d Llil d L ({27}, AI = 3/2); (5) 
5 = s Ll ,\ a d L { J2 QRl^qR) ({8}, AI = 1/2); (6) 

q=u,d,s 

6 = s L lM ]T g fl7At g fl ) ({8}, A/ = 1/2). (7) 

This set is sufficient for calculation of the CP-even parts of the amplitudes 
under consideration. To calculate the CP-odd parts , it is necessary to add 
the so-called electroweak contributions originated by the operators Oj, 0$: 

Or = \s llx {\ + 75M E - 75)?) (AI = 1/2, 3/2); (8) 

q=u,d,s 

O s = -12 £ e q (s L q R )(q R d L ), e q = (H -I, -I) (A/ = 1/2, 3/2). (9) 

q=u,d,s 

The coefficients 05-3 have the imaginary parts necessary for CP violation. 

In the case of non-linear realization of chiral symmetry, the bosonization 
of these operators can be done using the relations [4] 

%(l+n)qk = -^F w r(u-± i &U) , (10) 



fc7„(l + 7 5 )? fe = *[( W C/t - U (d^) - (m(9,C/t) _ (d,U)m)] kj 

(11) 

where 

-FV ^ 93 MeV, A 1 GeV, r = 2m^/(m u +m d ), m = diag(m u , m d , m s ). 

(12) 
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where a 3 is an arbitrary number and 
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(13) 



The PCAC condition demands = [5] and we adopt this condition, bear- 
ing in mind that, on mass shell, the values of the mesonic amplitudes are 
independent of a 3 . 

Using also the relations between matrices in the color space 

<W = \8ft} + \\ a s \} 

\a \7 _ 16 XaSTI 1 \a \7 
A f3 A 6 - ~9°S°I3 ~ 3 A <5 A /3 

and the Fierz transformation relation 

S7 M (1 + 7 5 )rf • ^(1 - 7 5 )g = -2s(l - 7 5 )g • q(l + 
and representing MiK — > 27r) in the form 



M{K\ -> 7T+7T-) = A e i5 ° - A 2 e i<52 , 
M(X 1 ° -> 7T°7r ) = A e l<5 ° + 2A 2 e i<52 , 

M(X + -> tt+tt ) = --A 2 e l52 , 
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(14) 
(15) 

(16) 



where <5 and 5 2 are the 5 -wave shifts of nn scattering in isotopic spin 7 = 0,2 
channels, we obtain 



A = GpF n sin 9 C cos 9 C 



rn 



K 



mlr 32 . , ( 

^—[ci - c 2 - c 3 + y/3(Rec 5 + ilmc 5 )]; (17) 



A 2 = Gf-^tt sin 6q cos 6*c 
where 



(18) 



C5 = c 5 + — c 6 , c 7 = c 7 + 3c 8 . 
16 







2mi 



A 2 (m u + m d ) 2 ' 



(19) 
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The contributions from Cj0 7 into ReA and Imv4 are small because c 7 /c 5 
is proportional to the electromagnetic constant a and we neglected these 
corrections. From data on widths of K — > 2rr decays we obtain 

32 

c 4 = 0.328; ci - c 2 - c 3 + — /3Rec 5 = -10.13. (20) 

9 

At ci — c 2 — c 3 = —2.89 [3], [6] we obtain 

32 

— /3Rec 5 = -7.24. (21) 
9 

From the expression for A 2 , it is seen that the contribution of the operators 
7)8 is enlarged by the factor A 2 /m 2 K in comparison with the other operator 
contribution. The reason is discussed in Appendix. 
Using the general relation 



£ t = ie i(S2-s Q ) 



ImAn ImA? 



^0 



ReA ReA 2 . 

and the experimental value e' = (3.78 ±0.38) x 10~ 6 , we come to the relation 

_ Imcs / _ 24.36^") = (1.63 ± 0.16) x 10" 4 , (23) 

Rec 5 V ,,n Imc 5 / V ' V ' 

where Q V:V > takes into account the effects of K° — > -K Q rj{rj') — > Tr°Tr° transi- 
tions. 

Introducing the notation 

Imc5 ImAf 24.36 Imc7 

— = x ■ = — y f 24) 

Recs si ' 1 — Imcs 

and using 

(ImA t )M « S2S3 sin 5 = ^^f^ [7] (25) 
we can write Eq.(23) for Q VtV ' = 0.25 ± 0.08 in the form 

x(l-y) = 0.40 x (1± 0.22). (26) 

In the last two equations and 5 are the parameters of CKM matrix. The 
Eq.(26) depends on the variables x and y representing the contribution of 
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QCD penguin and relative contribution of EWP, respectively. To move far- 
ther, we are enforced to apply to existing theoretical estimates of one of these 
variables. 

In terms of notations in [8-10] 

= 5^/110(1-^). (27) 

According to [8] y « 0.3 and hence x = 0.57±0.12. But e'/e = 2.2- 10~ 3 
or by 30% is larger than the experimental value . 

In [10], the central value of y is y ~ 0.5 and, consequently, x = 0.80±0.18. 
This result looks as the reliable one. A very close result x = 0.71 ± 0.27 can 
be derived from the result (e'/e^jjyyyp = (—12 ± 3) x 10~ 4 [11] comparing 
it with the experimental value [e'/e)exp = (16.6 ± 1.6) x 10~ 4 . But it should 
be noted that the previous estimates of x were rather different. In particular, 
according to [12] x — 1.4 ± 0.28. An estimate of x can be extracted also 
from the papers [13-15] operating with different set of 4-quark operators Qi, 
where the combination CqQq corresponds to our C5O5. From the general 
representation 

C 6 (n) = z 6 (/ji) + (si + s 2 s 3 —?- cos 5) ■ y 6 (//) - is 2 s 3 -^- sin 5 ■ y 6 (//) 

C1C3 C1C3 

and the calculated magnitudes of y 6 and z 6 we find for x ~ ye/z 6 : 

x ^2 at A^ CD = 0.35 GeV, = 0.8 GeV, m t = 176 GeV [13]; (28) 

x = 2.8 at A Ms = 0.3 GeV, fi = 1 GeV, m t = 130 GeV [14]; (29) 

x = 5.5 at Aq CD = 0.3 GeV, fi = 1 GeV, m t = 170 GeV [15]. (30) 

Such difference of the theoretical estimates of x makes very desirable an 
investigation of CP-effects in K ± — > 7r ± 7r ± 7r T decays, where, contrary to 
K L — > 2n decays, the EWP contributions increase CP effects . 

3. Decay K± 7r ± 7r ± 7r T 

To leading p 2 approximation 
M(K + -> 7T + (p 1 )7r + (p 2 )7r~(p 3 )) = «[1 + + ^V(l +^ M ) + ...], (31) 
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where 



d-KM 



k = Gp sin 9 C cos 6cm 2 K c (3V2) 1 , 

32 ^ _ /3A 2 „\ 

— plmc5 + 4plmc7 
9 



2m 2 



A" 



/c , 



'KM 



32 
IT 



/9Imc5 + 8/9Imc7 



/(c + 9c 4 ). 



(32) 
(33) 

(34) 



The last two quantities represent the imaginary parts produced by the Kobayashi- 
Maskawa phase 5. 



\9 = 



Sw 2 

^(l + 9c 4 /c ), 



2m 



Y = (s 3 - s ) /ml, 



K 



Cq — C\ — C2 — C 3 



32 

c 4 + — /5Rec 5 
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■10.46. 



(35) 



(36) 



As the field is the complex one and its phase is arbitrary, we can 
replace K + by K+(l + ia KM ){^l + a KM )~\ Then 



M (K + -> 7T + (p 1 )7r + (p 2 )7r-(p 3 )) = k[1 + -<?y(l + i(&^M - %m)) + ...]. (37) 

Though this expression contains the imaginary CP-odd part, it does not lead 
to observable CP effects. Such effects arise due to interference between CP- 
odd imaginary part and the CP-even imaginary part produced by rescattering 
of the final pions. Then 

M(K + -> 7r + 7r + 7r) = «[1 + ia + \gY{\ + ib + i(b KM - a KM ) + ■■■} (38) 

where a and b are corresponding CP-even imaginary parts of the amplitude. 
These parts can be estimated to leading approximation in momenta calcu- 
lating the imaginary part of the two-pion loop diagrams with 

M(7r+(r 2 )7r-(r 3 ) -> 7r+(p 2 )7r-(p 3 )) = F~ 2 [(p 2 + p 3 ) 2 + (r 2 - p 2 f - 2m 2 ], 
M(7r°(r 2 )7T (r 3 ) - 7r+(p 2 )7r-(p 3 )) = F" 2 [(p 2 + p 3 ) 2 - ml], 
M( 7 r+(r 1 )7r+(r 2 ) - 7r+( Pl )7r+(p 2 )) = F" 2 [(n - Pl ) 2 + (n -p 2 ) 2 - 2m 2 ]. 



Then we find: 



a = 0.12065, 6 = 0.714. 



(39) 
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Using the definition 



\M{K± - 7r ± (p 1 )7r ± (p 2 )^(p 3 ))\ 2 ~ [1 + g^Y + ...] 
and the results of our calculation 

\M{K± - 7i ± ( Pl )7i ± (p 2 )^(p 3 )\ 2 ~ (1 + ab ± a(6 XM - a™))+...] 

1 + cr 

(40) 

we find 

P = # + ~ 9' _ ajbKM ~ Okm) {a , v 

At the fixed above numerical values of the parameters and Q ri ,-q' — 0.25 we 
obtain to leading p 2 approximation 

(R g ) 2 = 0.030^(1-14.9^) = = -(2.44±0.44)xl0- 5 a ; (l - 0A3±0m 
V 9jp Rec 5 V Im5 5 ; v J \ x 

(42) 



4. The role of p A and other corrections 

The corrections to the result obtained in the conventional chiral theory up 
to leading p 2 approximation are of two kinds. The first kind corrections are 
connected with a necessity to get explanation of the observed enlargement of 
S-wave I = 7T7T amplitude. The corrections of the second kind are the p 4 
corrections. As it was argued in [16], [17], the corrections of both kinds be 
properly estimated in the framework of special linear U(3)l® U(3)r a model 
with broken chiral symmetry. The above mentioned enlargement of S wave 
in this model is originated by mixing between the qq state and the gluonic 
state (G^,) 2 states. In such a model 

U = a + in 

where a is 3 x 3 matrix of scalar partners of the mesons of pseudoscalar 
nonet. The relations between diquark combinations and spinless fields are 
as given by eqs.(10) ,(11), but without the terms proportional to A~ 2 . Such 
contributions in a model appear from an expansion of the propagators of 
the intermediate scalar mesons. The parameter A 2 in this model is equal to 
difference ^a (980) — m n- The strength of mixing between the isosinglet a 
meson and corresponding gluonic state is characterized by the parameter £. 
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If the p 2 approximation gives 

(k) p 2 = 1.495 x 10" 6 , (g) p 2 = -0.172 (43) 

instead of 

(«)exp = 1-92 x 10" 6 , (<?) eX p = -0.2154 ± 0.0035, (44) 

the corrected values of these CP-even parameters of K + — > 7r + 7r + 7r~ ampli- 
tude are closer or practically equal to the experimental ones [16]: 

(/t) ( p2 +p 4. ?= _ .22 5 ) = 1-73 x 10 ~ 6 , (50(p2+p4. 5= _ o .225) = -0.21. (45) 

More information on the parameter £ can be found in [16], [17]. The expres- 
sions for the corrected ttti — > tttt amplitudes are presented in [17]. 

Calculating the CP-even imaginary part of the K ± — > 7r ± 7r ± 7r T amplitude 
originated by two-pion intermediate states, we obtain 

a(p 2 + p 4 ; £ = -0.225) = 0.16265, (46) 

b(p 2 +p 4 ; £ = -0.225) = 0.762. (47) 

An estimate of the parameter a can be obtained also without any cal- 
culations using the experimental data on the phase shifts of nn scattering 
8q, 8q, 8\. According to definition (38) a is a phase at s 3 = s . The mean 
value of the squared energy of ty + ti~ system is 

1[(P1+P3) 2 + (P2+PS) 2 }=S +-^^. 

Consequently, a is a phase shift of 7r + 7r~ scattering at y/s = But 
the only significant phase shift at y/s = y/s^ is 8®. The rest phase shifts 
are very small: |5q( s o)| < 1.8° and 8l(s ) < 0.3° [18]. Then, according to 
Eq.(38), a « tan5°(s ), or a = 0.13±0.05, if 5g(s ) = (7.50±2.85)° [19] and 
a = 0.148 ± 0.018, if 8$(s ) = (8.4 ± 1.0)° [20]. These results coincide inside 
the error bars with the result (46). The corrected magnitude of R g is 

W; ^-o, 25 ) = 0.039^J (l - 11.95gg) = 

= -(3.0 ± 0.5) x I0~ 5 x (l - Miyogs ) , 1 ; 

This result is by 23% larger in absolute magnitude than that calculated in the 
leading approximation. Therefore, we come to conclusion that the corrections 
to the result obtained in the framework of conventional chiral theory to the 
leading approximation are not negligible (23%), but not so large, as it was 
declared in [21]. 
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5. Conclusion 



From Eqs.(22), (26) and (48), it follows that EWP contributions diminish 
e'/e and increase R g . The EWP corrections cancell one half of the QCD 
penguin contribution into e'/e at x — 0.8 and cancell 80% of QCD penguin 
contribution at x — 2. In both cases e'/e is the same. 

In the case of K ± — > 3n decays, the direct influence of EWP corrections 
themselves on CP effects is not so crucial as in Kl — > 2n decays. But if a 
cancellation between the contribution of QCD and electroweak penguins in 
e'/e is large, the factor x in Eq.(48) is also larger than 1. So, for x = 2 , the 
predicted R g must be 2.5 times larger than at x — 0.8. 

Therefore, measuring R g , one obtains a possibility to determine the true 
relation between QCD and EWP contributions into CP violation in kaon 
decays. 



Appendix 



Here we explain why an expansion of the amplitudes originated by elec- 
troweak penguin diagrams begins from the term, independent of momenta 
and masses of the pseudoscalar mesons. 

The operators 7j8 can be expressed in terms of colorless diquark combi- 
nations in the form 



7 = -s(l - 7 5 )m • w(l + l 5 )d - -0 5 , 
Using eq.(10), we find 

F 2 r 2 



8 = 30 7 . 



On 



-U 2 iU* 3 + (terms proportional to ■p i (mj). 



(A.l) 



(A.2) 



Omitting the terms proportional to derivatives of U and taking in Eq.(12) 
a 3 = 0, we find for the parity-even transitions 



(07 



\P— even 



Fir 



X 



3k + %) + « + K° +■■■}■ 



(A.3) 



This expression does not contain the direct contribution to K + — > 3n decays, 
but thanks to the term n~K + , the independent of p| (m|) part of K + — > 3n 



10 



amplitude arises. In p 2 approximation 

< 7I+( Pl )7I+(p 2 )7I-( P3 )\0 7 \K+(k) > = 



Fir 



si+s 2 -2m|. _ s 1 +s 2 -ml-m 2 K 



F 2 



F 2 



(A.4) 



2 ' 



where Sj = (k — pi) 2 and the first term in the brackets describes the n + (k) — > 
7r+ (Pi) 7r+ (j'2) 7r ~(P3) transition. The second term describes the transitions 
K + (k) ^ K + (p h2 )ir + (p 2A )ir-(p 3 ). 

Therefore, the operator O-j violates the rule, according to which an ex- 
pansion of the mesonic amplitudes begins from the terms proportional to 
p 2 {m 2 ). 

It may seem to one that removing the non-diagonal term — y7r~i^ + from 
the effective Lagrangian by redefinition of K + and tt~ fields [22], the problem 
with the constant contribution could be solved. But this is not so. 

In our case, the mass part of the effective Lagrangian contains, in partic- 
ular, the combination 



m 2 K K + K- 



P2 2 



(A.5) 



where 7 = V2G-p sin^Q cos 9qCj. The transformations 



7T- -> 7i- + PK~, K+^K + -(3n+, 

7T+ -> 7T+ + K- ~^K~ - (3*11- 



(A.6) 



with 



/3 = 7*F>72(m; 
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(A.7) 



remove the non-diagonal terms in the linear in 7 approximation. But the 
effective Lagrangian of strong interaction generates the sum of the amplitudes 

< 7l + ( Pl )7l + (p 2 )7l-(p 3 )\7V + (k) > + < K + ( Pl )7r + (p 2 )7T-(p 3 )\K + (k) > + 

<K+(p 2 )n+( Pl )n-(p 3 )\K + (k) > 

(A.8) 

which after the transformation (A. 6) generates the amplitude 



< n+( Pl )n+(p 2 )n-(p 3 )\0 7 \K + (k) >-- 



7* 



si+S2-2m* 
F 2 



s\+S2-m 2 T -m 2 K 



F 2 



2 ■ 



(A.9) 
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We have reproduced the result (A. 4). The contribution of the operators O-j^ 
to the leading approximation does not depend on pj{mf). 
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